Nonclassical correlation in NMR quadrupolar systems 
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The existence of quantum correlation (as revealed by quantum discord), other than entanglement 
and its role in quantum-information processing (QIP), is a current subject for discussion. In partic- 
ular, it has been suggested that this nonclassical correlation may provide computational speedup for 
some quantum algorithms. In this regard, bulk nuclear magnetic resonance (NMR) has been suc- 
cessfully used as a test bench for many QIP implementations, although it has also been continuously 
criticized for not presenting entanglement in most of the systems used so far. In this paper, we report 
a theoretical and experimental study on the dynamics of quantum and classical correlations in an 
NMR quadrupolar system. We present a method for computing the correlations from experimental 
NMR deviation-density matrices and show that, given the action of the nuclear-spin environment, 
the relaxation produces a monotonic time decay in the correlations. Although the experimental 
realizations were performed in a specific quadrupolar system, the main results presented here can 
be applied to whichever system uses a deviation-density matrix formalism. 
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I. INTRODUCTION 

Since the birth of quantum information theory (QIT), 
entanglement has been considered as a key resource 
for the processing of information at a quantum level. 
However, it is known that quantum correlation may be 
present even in separable states. This nonclassical cor- 
relation can be quantified, for example, by the so-called 
quantum discord pQ, which is a measure of a "gap" be- 
tween quantum- and classical- information theory. Some 
other measures for both classical and quantum correla- 
tions contained in a bipartite quantum state were pro- 
posed in literature Nonclassical correlation may 
exist in almost all quantum states [10) . and it has been 
theoretically [TT] and experimentally [T3] demonstrated 
that it may provide computational speedup in a model of 
quantum computation. Also, it was suggested [13] that 
the speedup of some quantum algorithms may be due to 
quantum correlation of separable states. 

The study of the nonclassical aspects of a correlated 
quantum system, especially the aspect revealed by quan- 
tum discord, received a great deal of attention in re- 
cent scientific literature [HI Q3H2Z]. In particular, the 
action of decoherence on this correlation was studied 
by taking different kinds of environmental interactions 
[151 , 123H251 128] into account. Until now, only two exper- 
imental measurements of such a nonclassical correlation 
have been performed [TH [53] . In Ref . [T2] , by means 
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of an optical architecture, the authors implemented the 
so-called deterministic quantum computation with one 
qubit [30J for the trace estimation of an unitary matrix. 
In such a non-universal quantum-information-processing 
(QIP) model, entanglement is not a necessary resource 
for obtaining a computational speedup (in comparison 
with the best classical algorithm). The authors of Ref. 
[29] demonstrated, also using an optical setup, the sud- 
den change in the decay rates of classical and quantum 
correlations, theoretically predicted in Ref. [T5] . 

Nuclear magnetic resonance (NMR) systems have been 
extensively used as a method for implementing and 
benching test QIP ideas (3JJ [32]. The main feature of 
the technique is the excellent control of unitary trans- 
formations provided by the use of radio-frequency (rf) 
pulses, which result in unique methods for quantum state 
manipulation [33 , 34J and generation of protocols for pro- 
cessing quantum information |31 [ I35 [ [36] . For example, 
algorithms such as Deutsch-Jozsa, Shor, and Grover were 
successfully tested using NMR systems [3TJ [35]. How- 
ever, most of these achievements were performed in bulk 
samples (bulk NMR) by using the so-called pseudop- 
ure states, for which the existence of entanglement was 
promptly ruled out [37H39] , This last fact has led to the 
questioning of the quantum nature of NMR implementa- 
tions for QIP. On the other hand, as suggested in Refs. 
[T3J[37], the existence of quantum correlation, other than 
entanglement in NMR systems, may be the reason for 
the success of most bulk NMR implementations. More- 
over, as mentioned before, quantum correlation of sepa- 
rable states may provide computational speedup in some 
tasks. 

To our knowledge, here, we present the first demon- 
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stration that both quantum and classical correlations in 
NMR systems can be determined from experimentally 
detectable deviation-density matrices. We use an NMR 
quadrupolar system to theoretically and experimentally 
study the existence of quantum and classical correlations 
in a two-logical-qubit composite system. Our results 
show that such a nonclassical correlation can easily be 
created and can be manipulated at room temperature. 
We also investigate how these correlations are degraded 
by the effect of the environment on the nuclear spin. In 
our system, the quantum aspect of correlation decays 
monotonically within the relaxation time of the deviation 
matrix. Our experimental implementations were carried 
out by using an effective two-qubit representation of an 
23 Na nuclear spin (I = 3/2) in a lyotropic liquid-crystal 
[sodium dodecyl sulfate (SDS)] sample at 26° C. The 
achievements presented here, together with the discus- 
sion about the role of quantum correlation (of separable 
states) in QIP, reinforce the relevance of the NMR tools 
in this scenario. 

This paper is organized as follows. In Sec. II, we briefly 
review the description of the NMR quadrupolar system. 
Section III is devoted to presenting the method for deter- 
mining the quantum and classical correlations from the 
experimentally accessible NMR deviation-density matrix. 
Also in this section, we present a theoretical model to 
describe the relaxation dynamics in our system. The ex- 
perimental results are presented in Sec. IV; and, in Sec. 
V, we present a summary and a discussion. 



II. DESCRIPTION OF THE SYSTEM 

In NMR systems, at room temperature, the energy 
gaps between the levels of the system SE — huj L are much 
smaller than the thermal energy e = /kj£,/2/cbT ~ 10~ 5 
(we added, for convenience, a 1/2 factor in the definition 
of e). This implies that a typical NMR-system density 
matrix can be written, in the high-temperature expansion 

eh ma, as 



the following way: 
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where 1 is the 2" x 2™ identity matrix, with n being the 
number of effective logical qubits in the QIP terminology, 
Z w 2" is the system-partition function [31 j . and Ap is 
the traceless deviation-density matrix. In general, the 
state described in Eq. |l| is a mixed state that does 
not possess entanglement 3H1 13S] ■ However, as we will 
show in what follows, it might have quantum correlation 
(of separable states) that can be used in QIP protocols 

num. 

Any manipulation — such as state preparation, state 
tomography, qubit rotations, and so on — is performed 
only onto the deviation-matrix Ap. For example, a se- 
quence of rf pulses, which can be represented by a uni- 
tary transformation U, changes the density operator in 



UpUl = —l + eUAprf. 



(2) 



By a suitable adjustment of each rf pulse length, phase, 
and amplitude, very fine control over the density matrix 
populations and coherences (diagonal and off-diagonal el- 
ements, respectively) can be achieved. Together with 
proper temporal or spatial averaging procedures [41] and 
evolution under spin interactions, the rf pulse can be spe- 
cially designed to prepare the system in all two-qubit 
computational base states as well as its superpositions, 
starting from the thermal-equilibrium state (Boltzmann 
distribution) [42j [43] . The effect of the environment on 
the spin system is to induce relaxation in such a way that, 
after characteristic times, the populations return to the 
equilibrium distribution and the coherences vanish. 

The purpose of this paper is to theoretically and exper- 
imentally investigate the presence of quantum and classi- 
cal correlations and their degradation due to decoherence 
for different initial states of a two-qubit representation of 
the nuclear spin system. To achieve this goal, we use a 
spin 1 = 3/2 quadrupolar NMR system. In the presence 
of a static magnetic field, due to the Zeeman splitting, 
nuclei with / > 1/2 can be described by a (21 + l)-level 
system with equally spaced energy levels. However, such 
nuclei also possess quadrupole moments, which interact 
with the electric-field gradient (EFG) produced by the 
charge distribution in their surroundings, the so-called 
quadrupolar interaction [31 j . When the Zeeman inter- 
action is much stronger than the quadrupolar one, the 
latter can be treated in the framework of first-order per- 
turbation theory, and the system's Hamiltonian turns out 
to be 



H = -Hu) L I z + hbJr 



{Ml 



(3) 



where u)l and ujq are the Larmor and the quadrupole 
frequencies, respectively (\u>l \ I w q|)- The spin-nuclear 
operator is characterized by its z component I z and its 
square modulus I 2 . The first term of Eq. ^ describes 
the Zeeman interaction while the second one accounts for 
the static first-order quadrupolar interaction [4"4"l - |4ll)] . 

This system can be regarded as an effective two-qubit 
quantum system that may be used for QIP [32 1321 13S1 
Fffl FIT] . Figure 1 shows a sketch of the four energy levels 
in the case I = 3/2, characterized by the z component 
of the nuclear-spin magnetization whose projections are 
m = -3/2, —1/2, 1/2, 3/2. Figure 1 also displays the 
logical indexing of the angular momentum basis [i.e., the 
computational basis (for the effective two-qubit system)] 
as well as the three single quantum transitions that can 
be directly detected after applying a non selective n / 2 rf 
pulse onto the system (equilibrium spectrum) . Although 
the implementation of the qubit in a quadrupolar system, 
as considered here, is not so obvious, many experiments 
proved that logic gates as well as quantum algorithms 
can be realized in such an effective system [JHl EH] • 
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FIG. 1: Sketch of the four energy levels in the spin / = 3/2 
system (characterized by the z component of the nuclear- 
spin magnetization) for the Zeeman and the quadrupolar cou- 
plings. The indexation of the effective two-qubit computa- 
tional basis is also displayed. 



III. QUANTUM AND CLASSICAL 
CORRELATIONS AND THEIR DYNAMICS 
UNDER DECOHERENCE 

Let us now briefly review two proposed measures to 
quantify the quantum and classical aspects of the cor- 
relation contained in a bipartite quantum system. The 
nonclassical correlation in a composed state can be quan- 
tified by quantum discord [T], and its classical counter- 
part is given by the Henderson- Vedral measure [5] . In 
what follows, we will also discuss a measure of quantum 
correlation that can be regarded as a symmetrical ver- 
sion of quantum discord and its classical coun- 
terpart [9] . A suitable expansion of the correlations 
measures, in terms of the experimental measurable de- 
viation matrix [cf. Eq. will then be obtained. Fi- 
nally, we present the dynamic equations governing the 
evolution of the density operator of a spin 1 — 5/2 sys- 
tem under the action of a typical spin environment in an 
NMR quadrupolar scenario, using the Redfield formalism 
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Measures of correlations 



The mutual information is a measure of the correla- 
tion between two random variables A and B in classical 
information theory [52] 



I C (A:B) = H{A) + H{B) - H(A, B) 



(4) 



where H(X) = — J2 X Px log 2 Px is the Shannon entropy of 
the variable X with p x being the probability of variable 
X assuming the value x. By means of the Bayes rule, we 
can rewrite Eq. Q in an equivalent form 



J C (A:B) = H(A) - H(A\B), 



(5) 



with H(A\B) = - J2 a ,bPa, b log 2 p a |b (p a \b = Pa,b/pb) be- 
ing the conditional entropy, which represents the lack of 
knowledge of variable A when we know (by means of a 



measurement) the value of variable B. It is clear that, in 
the classical case, we always have I C (A:B) — J C (A:B) = 0. 

In the realm of QIT, while the extension of Eq. Q for 
a bipartite quantum state {pab) is trivially obtained as 



A:B) 



S(p A ) + S(p B ) - S{pab), 



(6) 



the extension of Eq. ([5| is not so straightforward. Here 
S(p) = — Tr( ( olog 2 p) is the von Neumann entropy, and 
Pa(b) = r ^ T B(A)(pAB) is the reduced-density operator of 
partition A(B). The quantum mutual information given 
in Eq. ^ is a measure of the total correlation (includ- 
ing the classical and the quantum ones) contained in a 
bipartite quantum system [5, 56]. The quantum exten- 
sions of Eqs. Q and ^ are not equivalent, in general; 
and this inequivalence relies on the distinct nature of 
both quantum and classical measurements. While the 
classical measurement can be chosen such that it does 
not disturb the system to be measured; in the quantum 
case, the measurement process may affect the system in a 
fundamental manner. This observation leads Ollivier and 
Zurek pQ to propose a measure of quantum correlation, 
named quantum discord, given by the difference 
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A:B ) 



max J{pa-.b) 



(7) 



where J{pa-.b) = S(p A ) - Sr n? j (pa\b), wi th 

S{n B } (pa\b) = (pa) being a quantum exten- 

sion of the classical conditional entropy T-L{A\B). The 
reduced state of partition A (p^), after the measure- 
ment is performed on partition B, is given by 

Pa = Tr s {(lA®nf)p AB (l A <SUf)}/qj, with q 3 = 
Ttab [(iA®nf) Pab\ (1a is the identity operator for 
partition A). The quantum discord Eq. ^ is computed 
by an extremization procedure over all possible complete 
sets of projective measurements {n^} over subsystem B. 

Due to the distinct nature of quantum and classical 
correlations, we can assume that both correlations add 
up in a simple way to give the quantum mutual informa- 
tion. Therefore, the classical counterpart of the quantum 
discord may be defined simply as [13 125] 



C(pab) = Z(pa-.b) ~ T^(pab), 



(8) 



that is identical to the Henderson- Vedral definition of 
classical correlation [21 [3] . 

In some circumstances the quantum discord, Eq. ([7|, 
and also its classical counterpart, Eq. ([8]), may be asym- 
metric with respect to the choice of the partition to be 
measured (see Refs. [25], 126] for a related discussion con- 
cerning symmetric measures of correlations). Since the 
states considered in this paper are, in general, asymmet- 
ric by an interchange between the two partitions we will 
use symmetrized versions for the measures of both quan- 
tum and classical correlations. Alternatively, the clas- 
sical correlation in a composite bipartite system can be 
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expressed as the maximum classical mutual information 
that is obtained by local measurements over both parti- 
tions of a composite state [HI [5] 



IC(pab) = max [I C (A:B)} , 
{nf®nf } 



(9) 



where I C (A:B) is the classical mutual information defined 
in Eq. Q for the probability distributions that result 
from the quantum-measurement process. The extrem- 
ization in Eq. ([9| is taken over the set of local projective 
mcasur ements \U.f ^ n f } 

over both subsystems. Since 
the quantum mutual information quantifies the total cor- 
relation, a symmetrized measure of the quantum corre- 
lation can be defined as E51 |2"61 



Q(pab) = Z(pa-.b) - IC(pab)- 



(10) 



For composite states of two qubits with maximally mixed 
marginals the quantum discord Eq. Q is identical to its 
symmetrized version Eq. (10), i.e., T>(pab) — Q(pab) 
[and also IC(pab) = C(pab)] 
in general [25J [2~ 



However, it is not true 
We can regard Eq. 



(10) as a sym- 



metrical version of the quantum discord. Indeed such a 
symmetrical quantifier also reveals a departure between 
the quantum and the classical versions of information 
theory. 



Correlations and the deviation matrix 



By replacing Eqs. ( 12 ) and ( 13 1 into Eq. ([6| and by 
disregarding terms containing high powers in e, we obtain 



In 2 



(2 Tr Ap 2 - Tr Ap 2 A - Tr Apj 



(14) 



To compute the classical correlation we must obtain 
the measured density operator, which is given by 



V = J2 ( n f ® n f ) p ( n f ® u j 



i 



eAij, 



(15) 



where we defined the measured deviation matrix as 
Ar) = ( u ? ® n f ) A P ( n f ® nf ). For a two-qubit 
system, the complete set of projective measurements is 
given by {IT) = |6j) (0*| , j =||, JL , ft = A,B}, where 

OfA = coefflfcllOk) + e^ fc sin(0 fc )|l fc ) and 



lei 



) 



= cos(0 fc ) |0 fc ) + e^ fc sin( 

e-*** sin(0 fe ) |0fe) - cos(0 fe )|lfc) with < 6 k < tt and 
< < 2-7T. {|0fc) , |lfc)} is the computational basis of 
the logical qubit of partition k. We note that the corre- 
lations quantifiers, presented here, have the same values 
for simultaneous or successive measurements performed 
on each partition. 



The same reasoning that results in Eq. ( 14 ) leads us to 



the following expression, in terms of the measured devia- 
tion matrix, for the mutual information of the measured 
state, 



For our purposes, the initial NMR density matrix can 
be written in terms of the deviation matrix as 



P = 4 + eA P- 



(11) 



The deviation matrix, Ap, can be experimentally recon- 
structed using a set of rf pulses and readouts [quantum- 
state tomography (QST)] [3TJ[32]. The parameter e may 
be estimated from the Zeeman and thermal-energy ra- 
tios (in our experiment, at room temperature, e <C 1). 
Since all correlations present in the state p comes from 
the deviation matrix Ap, it is convenient to express the 
measures of correlations as functions of Ap. To do this, 
we will expand the von Neumann's entropy in powers of 
the small parameter e as follows: 



In 2 



2Tr Ar] 2 - Tr (Ar/A) ~ ^{AtjeY 



and, thus, for the classical correlation, 



K{p) 



max /, 

{nf®nf } 



(V) 



(16) 



(17) 



where At]a(b) — T^BtA) Ar l 1S the reduced measured de- 
viation matrix of partition A(B). The quantum cor- 
relation in the composed two-qubit system can be di- 
rectly computed from Eqs. ( [l4| and (17) as Q(p) — 
l(p)-JC(p). 



Action of the environment on the deviation 
matrix 



S(p)=2h-^TrAp 2 j +••., (12) 

where we used the fact that Tr Ap = 0. The reduced- 
density operators reads Pa(b) = T^BfA) P = 1a(b)/2 + 
eApA(B), with Apa(b) — ^b(A) Ap. We observe that 
because Tr Ap = we have Tr Ap^ = Tr Ap B = 0. Thus 
the marginal entropies are straightforwardly obtained as 

2 

S(p A(B) )=l-^TrAp A{B) + .... (13) 



As commented on Sec. II, the quadrupolar coupling 
emerges from the interaction of the nuclear quadrupole 
moment with the EFG generated by the surrounding en- 
vironment. Internal molecular or atomic motions cause 
random fluctuations in the EFG, which introduces noise 
into the system, which leads to the relaxation process 
that causes decoherence and energy dissipation. The de- 
pendence of the random EFG fluctuations on the molec- 
ular motions is established by the spectral densities that 
encode the motion characteristics, such as geometry and 
correlation times [65] . 
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FIG. 2: (Color online) 23 Na NMR equilibrium spectrum of the 
oriented liquid-crystal SDS sample. The frequencies uoi, V12, 
and V23 correspond to the transitions 1 00} -o |01), |01) <-> |10), 
and 1 10} o 1 11), respectively. 
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FIG. 3: (Color online) Scheme of the pulse sequence used in 
the experimental procedure. The first step is the state prepa- 
ration through a sequence of SMPs [55]. Next, we leave the 
system to evolve only under the action of the decoherent en- 
vironment during a variable period of time t. The last step 
is the nonselective tomography pulse and the observation of 
the free induction decay (FID) signal, which enable us to re- 
construct the deviation-density matrix Ap through QST J5JJ . 
The 7T pulse was added in the middle of the free evolution 
period to refocus the Bo-field inhomogeneities. 



Our experimental implementations were carried-out 
using 23 Na nuclear spins (/ = 3/2) in a lyotropic liq- 
uid crystal at 26° C. In this system, the relaxation of the 
23 Na nuclear spin can be well described by considering a 
pure quadrupolar mechanism (the relaxation is produced 
exclusively by the EFG fluctuations), so it can be repre- 
sented by three reduced spectral densities Jo, Ji, and J 2 
at the Larmor frequency uil- The specific model that 
relates the molecular motions and the reduced spectral 
densities can be found in Ref. [55] ■ By applying the Red- 
field formalism [51, to our system, enables us to obtain 
the dynamic evolution of the deviation-density matrix 
elements Ap^- (i, j running from to 3) as [55] 



Apoi (*) = \ [Ap^ 



(ApS 



01 



+ a,4 

A/4) e" 



■2CJ 2 tl e -C(J + Ji)t 



(18a) 



Ap? 



13 



+ (Ap° 02 - Ap? 3 ) e' 2CJ ^] e 



-C(J +J 2 )t 



(18b) 



2 M 2 

(A/4 - 



+ Ap° 3 
Ap? 3 



e -2CJitl e -C(J + J 2 )t 



(18c) 



Ap 23 (t) 



~[Ap 01 + Ap° 23 



(Ap° 01 -Ap° 23 ) e 
°8a 



-2CJ 2 tl e -C(J + Ji)t 



Ap 03 (t) = Ap° 03 e- c ^+ J ^, 



A Pl2 (t) = Ap° 12 e~ c ^ +J ^ 
l r 
~ 4 



Apoo (t) 
Apn (t) = I 
Ap 22 (t) 

Ap 33 (t) 



Rx e" 



2C(Ji+J 2 )t 



R 2 e~ 2Cj2t -R 3 e - 2CJl *l 



l r 

4 



Ri e 



-2C(Ji+7 2 )t 



+i?2 e~ 2Cj2t -R 3 e~ 2CJlt ] 



1 r 

4 



i?i e 



-2C(Ji + J 2 )t 



i? 2 e~ 2Cj2t +R 3 e~ 2CJlt ] 



1 r 

4 



fix e - 2C (^+^)' 
+i? 2 e- 2Cj2t +R 3 e - 2CJlt ] . 



(18d) 
(18e) 
(18f) 

(18g) 
(18h) 
(18i) 
(18j) 



In Eqs. (18a|-(18j), the superscript refers to the 
initial value of each deviation-matrix element and Ri 
(i = 1, 2, 3) are constant coefficients. The labels of the 
deviation-matrix elements {0, 1, 2, 3} refer to the compu- 
tational basis of the effective two-qubit system ordered in 
the usual way {|00) , |01) , 1 10) , |11)} . The parameter C 
is proportional to the quadrupolar coupling ujq and can 
be obtained from the equilibrium NMR spectrum (dis- 
played in Fig. 2) [371 [55J. To fully describe the system 
relaxation, it is also necessary to determine the spectral 
densities Jj (i = 0, 1, 2) as well as the coefficients Rj 
(j = 1, 2, 3). This can be achieved by preparing an ini- 
tial state by using the technique of strongly modulated 
pulses (SMP) [59 61, together with temporal averaging, 
where all Ap coherences are non zero (the full superposi- 
tion state), by letting it evolves under the action of relax- 
ation during a time period r, and then by measuring each 
Ap element using QST 61 . By repeating this procedure 
for different values of r, the relaxation dynamics for each 
Ap clement is experimentally measured. In Fig. 3, we de- 
pict the pulse-sequence scheme of the experimental pro- 
cedure used. We note that Eqs. ( 18a[ )-( 18j I arc valid for 
whichever initial state of the system, with the asymptotic 
state being the equilibrium one, p = 1/4 + (21 z- 

Equations ( 18a )-( 18j ) can be combined to provide 
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single-exponential functions in the following way 



Ap 01 + Ap 23 = (Ap^ + Ap° 3 ) e -C(Jo+J 1 )t i (19a) 
Ap„ 2 Ap 13 = (Ap° 2 + A P %) e -c(Jo+J,)t ) (19b) 
Ap 12 = A^ 2 e" C(Jl+J2)t , (19c) 



Apoo 



Apn + Ap 22 — Ap 33 = Rie 



, (19d) 

t 


Re Ap 
0.6 r 

0.3 


' (19e) 


0.0 ■ 
-0.3 


• (19f) 


-0.6 
-0.9 ; 



For convenience, in Eqs. (19a)-(19f) we have suppressed 
the temporal dependence of the deviation-matrix ele- 
ments Apij(t). Thus, by fitting the experimental evolu- 
tion of the deviation matrix elements, it is possible to de- 
termine the mean values of the reduced spectral densities 
Ji by using Eqs. ( 19a )-( 19c ) and to determine the param- 
eters Rj by using Eqs. (19d)-(19f|. We also checked the 
consistency of the experimental data with the adopted 
relaxation model by comparing the values of J\ and J 2 
obtained from Eqs. 
Ref. 



(19a 



_ - fll9cD and Eqs. (|19d|)-(|l9f|) (see 
[58] for specific details). 
The relaxation dynamics of the system can also be 
described by the time dependence of the mean value of 
I z and I x>y [(I z }(t) - longitudinal relaxation, (I x ,y){t) - 



transverse relaxation] 
Eqs.(pa|-(pjl) as 



which can be calculated by using 



(I, 



{3Ap 00 - 



Apn - Ap 22 - 3Ap 33 } - 5 



2C 



- 2Cj2 * + i? 3 e- 2CJlt } 



(20) 



and 



(l X )(t)+i(ly)(t)=2A P21 

= Ap° ie - c ( Jl+j2 )*. (21) 



In Eq. (20), (I z )t represents (I z ) for the thermal- 
equilibrium state. The longitudinal relaxation is then 
characterized by the time constants tl 1 = (2CJi) _1 and 
tl 2 = (2C J2) 1 , and the transverse relaxation is charac- 
terized by the time constant r T = [C{ J\ + J 2 )]~ [rJ2"IlrJ3] . 

It is important to note that, in our effective two-qubit 
representation of the four-level system, the environment 
acts globally (i.e., it acts simultaneously on both logical 
qubits). Such global action of the environment can be 
observed in the form of Eqs. ( 18a| )-( [l8jj ), where we do not 
have distinct spectral densities for each qubit. We recall 
that we can consider this 7 = 3/2 system as an effective 
two-qubit one since we can manipulate the transitions 
between the energy levels of these logical qubits (i.e., the 
nuclear transitions of the quadrupolar spin) just as it is 
performed in the case of physical qubits |43| . 




FIG. 4: (Color online) Bar representation of the QST exper- 
imental data for the prepared deviation matrix: (a) real and 
(b) imaginary parts of the equivalent to the |X ran dom) pseu- 
dopure state; (c) real and (d) imaginary parts of the equiv- 
alent to the |$~) Bell-basis pseudopure state; (e) real and 
(f) imaginary parts of the equivalent to the |$ + ) Bell- basis 
pseudopure state. 



IV. EXPERIMENTAL RESULTS 

In our experiment, the sample of 23 Na nuclei in a ly- 
otropic liquid crystal, was prepared with 20.9 wt% of 
SDS (95% of purity), 3.7 wt% of decanol, and 75.4 wt% 
of deuterium oxide, by following the procedure described 
in Ref. [M]. The 23 Na NMR experiments were per- 
formed in a 9.4 T - VARIAN INOVA spectrometer us- 
ing a 7 mm solid-state NMR probe head at T = 26°C. 
Figure 2 shows the equilibrium spectrum for our sam- 
ple. From this spectrum we obtained the quadrupole 
frequency vq = 6wq/27t = 15kHz and the constant pa- 
rameter C = (12 ± 1) x 10 9 s~ 2 . By following the pro- 
cedure described in Sec. Ill C, we estimated the val- 
ues of the spectral densities as Jo = (17 ± 4) x 10 _9 s, 
Ji = (3.0 ±0.5) x 10~ 9 s, and J 2 = (3.4 ±0.5) x 10" 9 s. 
Consequently, the time constants that appear in Eqs. 
^ and (j^TJ are t l1 = (13 ± 4) ms, t L2 = (14 ± 4) ms, 
and tt = (13 ± 4) ms. Finally, the mean values obtained 
for the constant parameters Rj are R\ = —1.9 ± 0.2, 
R 2 = -7.9 ± 0.1, and R 3 = -4.5 ± 0.2. 

By means of numerically optimized rf pulses 



(SMP technique) 



and temporal averag- 
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FIG. 5: (Color online) Bar representation of the QST ex- 
perimental data for the prepared deviation matrix: (a) real 
and (b) imaginary parts of the equivalent to the Bell- 
basis pseudopure state; (c) real and (d) imaginary parts of 
the equivalent to the Bell-basis pseudopure state; (e) 

real and (f) imaginary parts of the equilibrium state. 
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FIG. 6: (Color online) Correlations in the effective two-qubit 
pseudopure states. The solid lines represent the theoretical 
model, and the marks are the experimental points. 



ing, we prepared the deviation-density matrices, 
which correspond to five distinct initial sates: 
an arbitrary random X-type pseudopure state 
(l^random)) and the four Bell-basis pseudopure states 

{|*±) = (|00)±|ll))/v^,|* ± >= (|01)±|10))/V2}. 
All of these initial pseudopure states have the following 
form, in the computational basis, for the corresponding 
deviation matrices, 1 



a / 
b e 
e* c 
f*00d 



(22) 



The five prepared initial-state deviation matrices are 
displayed as a bar representation in Figs. 4(a)-4(f) 
and Figs. 5(a)-5(d). The mean fidelity of the 
prepared states relative to the ideal ones is J- = 

1 /2 

(VPidcaiPproparodVPidoai) » 0.98. We note that, in- 
dependent of the initial state considered, the asymptotic 



1 In our experiment the null elements of the deviation-density ma- 
trix are zero within the experimental errors. 
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FIG. 7: (Color online) Correlations in the effective two-qubit 
pseudopure states. The solid lines represent the theoretical 
model, and the marks are the experimental points. 
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state is always given by the thermal equilibrium state 
whose deviation-matrix is Ap = 2I Z [displayed in Figs. 5 
(e) and 5(f)]. 

The reconstruction of the deviation-density matrix was 
performed through QST described in Ref. [61j. The 
relaxation dynamics of the elements was followed using 
the procedure described in Sec. Ill C. We computed the 
quantum correlation and its classical counterpart present 
in the effective two-qubit pseudopure state accordingly 
Sec. Ill B [i.e., they were computed through Eqs. (10 1, 

(fil 



and (17) by using the experimental reconstructed 
deviation-matrix for different evolution periods r = nAr 
(n = 1, ..,40)]. For each initial state, the free evolution 
depicted in the pulse sequence displayed in Fig. 3 was 
performed 40 times with an evolution-period increment 
of At = 1.5 ms, leading to a maximum evolution time of 
60 ms. The errors bars in Figs. 6(a)-6(c) and 7(a) and 
7(b) were estimated from the standard deviation of three 
or more realizations of the whole experimental procedure 
(sketched in Fig. 3) for each initial state (prepared via 
the SMP technique). 

We display the decoherence dynamics of correlations 
(classical, quantum and the mutual information) present 
in the two-qubit pseudopure states corresponding to 
Random), 1*+), I*"), and |$~) in Figs. 6(a)- 

6(c), 7(a), and 7(b), respectively. The theoretical curves 
obtained from the model presented in Sec. Ill C are 
in very good agreement with the experimental data. 
In these figures, we observe that the correlations de- 
cay monotonically (as an exponential law) under the ac- 
tion of the nuclear-spin environment. We also note that 
for the Bell-basis pseudopure states [|vl/ ± ) , l^), Figs. 
6(b) and 6(c), and 7(a) and 7(b)], the quantum corre- 
lation is typically greater than its classical counterpart 
Q(pab) > IC(pab)- The opposite situation occurs for 
the chosen |A ran dom) pseudopure state [Fig. 6(a)] where 
Q{pab) < K-(pab)- The former experimental observa- 
tion contradicts the early conjecture that the classical 
support of the correlation was believed to be equal to 
or greater than its quantumness for any quantum state 
[3] [31 |SJ [T7] . However, such an observation is consistent 
with previous theoretical predictions [Q3 [35] . It is also 
interesting to compare the dynamics of the quantum cor- 
relation with the relaxation of the system as described 
in Sec. Ill C. Note, for example, that the decoherence 
time of the quantum correlation for the state \X ran( i om ) 
is about 8 ms, the same order of magnitude for the re- 
laxation characteristic times tt, tl^, i~l 2 ~ 13 ms. The 
others initial pseudopure states exhibit a similar behav- 



fully implemented and have been tested in such a sce- 
nario [31] I37H39] . The debate about the quantum nature 
of NMR implementations for QIP [37] is renew due to the 
fact that separable states can exhibit nonclassical corre- 
lation, which may be responsible for the computational 
speedup in the NMR context [T3J . Therefore the investi- 
gations about such general correlation in an NMR system 
become quite relevant. 

In this paper we theoretically and experimentally stud- 
ied the dynamics of bipartite quantum and classical cor- 
relations in an NMR quadrupolar system under the ac- 
tion of the decoherence process (mainly caused by the 
EFG fluctuations). We reported the first observation of 
a nonclassical correlation in an effective two-qubit NMR 
system. We also provide an approach for computing 
the quantum and classical correlations in such a com- 
posed system from the experimentally accessible NMR 
deviation-density matrix. 

We have found that classical and quantum correla- 
tions decay monotonically in time, by following an ex- 
ponential law, which is in perfect agreement with the be- 
havior predicted by the theoretical model that was pre- 
sented. By depending on the initial state (random or 
Bell-basis pseudopure states), the relation between the 
amount of quantum and classical correlations changes. 
Our results show that the nonclassical correlation can be 
generated and can be manipulated in spin-pseudopure 
states of NMR systems at room temperature. In our 
experiments, the quantum aspect of correlation decays 
monotonically within the longitudinal and transverse re- 
laxation characteristic times. Although the experimen- 
tal realizations were performed in a specific quadrupolar 
NMR system, the methods developed here are quite gen- 
eral and can readily be applied to other NMR systems or 
whichever system uses a deviation-density matrix formal- 
ism. It is worth mentioning that, despite the fact that we 
prepared specific pseudopure states used in the context of 
the NMR QIP, the results are valid for any initial NMR 
state concerning its specific Ap. In other words, initial 
states obtained after simpler excitation schemes without 
temporal or spatial averaging may also have nonclassical 
correlations that show similar behavior. The results pre- 
sented here may shed new light on the quantum resources 
available in the NMR scenario for QIP implementations. 
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